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SYMBOLS 


B 


rigid body located at point P in which the vectors 1, j, k 
are fixed 


c 


( ) 


GJ 


cosine ( ) 

beam torsional rigidity, N-m^ 


H kinetic energy minus potential energy, N-m 

^ j ^ ]c unit vectors fixed in body B that, for the initial position of 

B, lie along the axes x,y,z and for the final position of 
B, the deformed state, lie along the axes X3, y^, 

K unit vectors along the reference axes x,y,z (fig. 1 ) 


L 


beam length, m 


L , distributed applied loads in the x,y,z directions, respec- 

tively, N/m 

M distributed applied moment about the x^ axis (i.e., an applied 

^ pitching moment) , N-m/m 

O(e^) a quantity that has order of magnitude less than or equal to 

third degree in the rotations 

p a generic point on the beam elastic axis 

r the distance along the deformed— beam elastic axis measured from 

the beam root, m 


) 

[T] 


sine ( ) 

the transformation matrix of direction cosines relating the 
axes X3, y^y z^ to x, y, z 

time, sec 


t^, t2 arbitrary times, sec 

U deflection of point P in the x direction, referred to as 

axial deflection, m 

Ug axial deflection of point P due to longitudinal strain of the 

elastic axis, a spatial quasi-coordinate defined as 

u + ( 1 / 2 ) f (v'2 + w'2)dx, m 

•^0 

V deflection of point P in the y direction, referred to as 

lead- lag deflection, m 


v 



w 


X, y, z 


^ 3 ’ ^ 3 ’ ^3 


^0 

6 

«<P 


deflection of point P in the z direction, referred to as 
flap deflection, m 

a set of Cartesian coordinates with origin at the beam root, the 
X axis being along the undeformed beam, m (fig, 1) 

a set of Cartesian coordinates with origin at point P with x 
remaining tangent to the elastic axis and and z^ along^ 

Principal axes for the cross section in the deformed state, m 

a rotation (flap) of the body B about the unit vector - j; the 
exact axis about which the rotation occurs is dependent on 
the sequence of angular rotations, rad 

angle with the x - y plane of the beam in the example 
problem, rad 

variational operator; also Dirac delta function 

virtual work per unit length done by nonconservative applied 
loads, N-m/m 

vector of virtual rotations of the reference frame X3, y^, z^ 
fixed in the beam cross section, rad 

a rotation (lead-lag) of the body B about the unit vector k; 
the exact axis about which the rotation occurs is dependent 
on the sequence of rotations, rad 

a rotation (pitch) of the body B about the unit vector i; the 
axis about which the rotation occurs is dependent on the 
sequence of rotations, rad 

pretwist angle, rad 

the third angle in the series lag-flap-pitch; it is the third 
angle required to rotate body B to the final position after 
first lag then flap rotations, rad 

the third angle in the series flap-lag-pitch; it is the third 
angle required to rotate body B to the final position after 
first flap then lag rotations, rad 

the k component of curvature of the beam in the example prob- 
lem, m”^ 

rotation per-unit-length vector of the reference frame x^, y , 
Z3; the 1 component is the torsion (angle of twist per unit 
length) and the j and k components are bending curvatures, 
m“^ 


vi 




( ). -V (vector) • i, j, k 

-*- > J > K. 

( j (vector) • I, J, K 
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SUMMARY 


Nonlinear expressions are developed to relate the orientation of the 
deformed-beam cross section, torsion, local components of bending curvature, 
angular velocity, and virtual rotation to deformation variables. These expres- 
sions are developed in an exact manner in terms of a quasi-coordinate in the 
space domain for the torsion variable. The entire formulation is independent 
of the sequence of the three rotations used to describe the orientation of the 
deformed-beam cross section. For more common cases in the literature in which 
one of the three rotation angles is used as the torsion variable, the result- 
ing equations depend on the choice of the three angles. Differences in the 
equations, however, are demonstrated to be in form only. The present deformed- 
beam kinematic quantities are proven to be equivalent to those derived from 
various rotation sequences by identifying appropriate changes of variable 
based on fundamental uniqueness properties of the deformed-beam geometry. 

This development helps to clarify the issues raised in the literature concern- 
ing the choice of the angles. The torsion variable used herein is also shown 
to be mathematically analogous to an axial deflection variable that has been 
commonly used in the literature. Both variables are quasi-coordinates in the 
space domain and have been used in derivations based on Hamilton’s principle, 
despite lack of rigorous justification. Rigorous applicability of Hamilton’s 
principle to systems described by a class of quasi-coordinates that includes 
these variables is formally established. 


1. INTRODUCTION 


1 . 1 Background 

To adequately model a helicopter blade, the deflections must be treated 
as moderately large and the equations of motion will be nonlinear. Such equa- 
tions are required for properly analyzing the stability and forced response of 
helicopter rotor blades in hover and in forward flight. To derive such equa- 
tions, it is necessary to first specify the geometry of the beam both in its 
undeformed state and in its deformed state at some particular instant in time. 
For typical beam theories, this involves expressing the position of a generic 


point P on the elastic axis, the orientation of a frame consisting of the 
axes normal to and along principal axes for the cross section at P, and any 
further deformations, such as warp of the cross section out of its nominal 
plane, to adequately specify the location of every material point in the beam 
(ref. 1). In the discussion of kinematics that follows we will focus on the 
variables used to describe the location of point P on the elastic axis and 
the orientation of the reference axes attached to the beam cross section at P, 
assuming that no deformation of the beam cross section occurs. This would, in 
general, imply a total of six variables; three deflection variables for the 
location of P and three angles^ for the orientation of the reference axes 
with origin at P with respect to a suitable space-fixed reference frame. 

For an Euler-Bernoulli beam, however, the cross section is assumed to remain 
normal to the beam elastic axis during deformation. Thus, two of the three 
angles can be eliminated by expressing them in terms of derivatives of the 
deflection variables. It is then necessary to express all the kinematical 
quantities in terms of the four remaining variables. 

It is clear that a material point can only occupy one position at a time. 
The beam geometry, by which we refer collectively to the position of the 
material points along the elastic axis and the orientation of the cross sec- 
tions in terms of direction cosines for the reference axes at P, is then 
unique at a specific time. When the beam geometry is specified, so is the 
strain, the torsion (or angle of twist per unit length), and the components of 
bending curvature along the principal axes of the cross section at P, since 
these also are geometric quantities. Although the values of these quantities 
are unique at a particular instant in time, the mathematical variables used in 
defining the geometry are not unique. In reference 2, Reissner shows, for 
example, that it is possible to formulate a large-deflection beam theory com- 
pletely independent of the choice of deformation variables. These equations 
are an extension of the Kirchhof f-Love equations (ref. 1) and illustrate the 
obvious fact that the geometric quantities describing the beam deformation 
under a given load, although the expressions may look different for different 
choices of variables, are determined by the loads alone. There are many possi- 
bilities, for example, in the choice of angles used to relate the orientation 
of the reference axes at P to a space-fixed Cartesian frame (fig. 1). In 
general, any change in orientation can be described by three successive rota- 
tions about three distinct, well-defined axes. These sequential rotations can 
either occur about space-fixed axes, or about axes that are fixed in an updated 
reference frame resulting from a previous rotation. The angles themselves are 
not unique, but values of the resulting direction cosines that define the 
orientation are unique. 


^It is possible to use variables other than angles to describe changes in 
orientation such as Rodriguez parameters (ref. 2). In this paper, however, 
we will restrict the discussion to the use of angles. 
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Figure 1.- Undeformed, pretwisted beam and reference Cartesian axes x, y, z. 


1.2 Discussion of Previous Work 

With the assumption that the beam cross section remains normal to the 
elastic axis during deformation, the exact expressions describing the deformed 
beam geometry have been developed and applied to a derivation of the equations 
of motion for a helicopter blade (refs. 3-5), In reference 3, the exact direc- 
tion cosines for the reference axes with origin at P were derived. These 
(ji]^ 0 ctlon cosines are invariant with the choice of angles used in describing 
the orientation. The direction cosines constitute elements of a transformation 
matrix [T] relating the reference axes with origin at P and a space-fixed 
set of Cartesian axes, the x axis of which is along the undeformed beam 
(fig. 1). 

In reference 4, the structural and inertial operators were derived, 
through second-degree nonlinearity, for the rotating beam. The bending curva- 
tures were developed in the appendix of reference 4, These expressions are 
exact and independent of the choice of angles used to describe the orientation. 
In reference 5, expressions for the aerodynamic loads are developed for the 
rotating beam. To derive expressions for the aerodynamic loads one must make 
use of components of each of the following vector quantities along the refer- 
ence axes with origin at P: the relative velocity of point P with respect 
to the fluid, the angular velocity of a rigid body B in which are fixed the 
reference axes with origin at P, and the virtual rotation of B (fig. 2). 

The first of these three items depends only on the [T] matrix, but the latter 
two must be developed separately. Only approximate forms of these quantities 
are given in reference 5 and the exact expressions, developed herein, have not 
previously been published. 

Torsion deformation in references 3-5 is described by use of a quasi- 
coordinate in the space domain, the elastic component of the angle of twist . 
The variable is referred to as a quasi— coordinate because it is related to 
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Figure 2.- Deformed-beam geometry. 


angles through integrals that cannot be evaluated in closed form. Hence, all 
three angles, no matter how they are chosen, can be eliminated from the deri- 
vation for an Euler-Bernoulli beam. Although the definition of the torsion 
variable is mathematically explicit in references 3-5, it is never stated 
therein that it is a quasi-coordinate. This may account for some of the dis- 
cussions that have appeared in the literature (refs. 6-8). 

To better understand these discussions it is helpful to consider a simple 
approximation. Sometimes an elastic beam is approximated by a system of 
spring-restrained rigid bodies, the simplest of which is a single rigid body, 
hinged at a point corresponding to the beam root to provide three angular 
degrees of freedom. Now, for a rigid— body system such as this, the arrangement 
of hinges presents the analyst with a choice of infinitely many distinct physi- 
cal spring arrangements — each of which constitutes a distinct, well-defined 
physical system. For example, in reference 3 differences in stability between 
two of these systems are identified. They each have only two degrees of free- 
dom: flap (out of the plane of rotation) and lead-lag (nominally in the plane 

of rotation). The only difference between the two systems is the hinge 
arrangement. By definition, lag-flap means that the flap hinge leads or lags 
with the blade; flap-lag means that the lead- lag hinge flaps with the blade. 
Thus, lead-lag motion is precisely inplane only for the lag-flap arrangement; 
hence, the word nominally is used in the above definition. For small angles, 
the equations are linear and the two systems have identical equations of 
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motion. Any degree of nonlinearity produced by moderately large angles will, 
however, produce differences in the blade pitch orientation for the two sys- 
tems for given flap and lead— lag deflections. These differences, although not 
large, produce different dynamic behavior and stability for the two systems. 

It is clear, however, that an elastic beam with a particular set of end con- 
straints is only one physical system and its kinematics, stability, and 
dynamic behavior are n^ subject to the analyst's choice of angles used in the 
derivation of equations that describe the system. While this may seem clear, 
there has been considerable discussion and some misunderstanding of this 
issue in the literature (refs. 6-8). 

The role of hinge arrangement in determining coupled flap-lag stability 
of a rotating, centrally hinged, rigid-beam model of a helicopter rotor blade 
in hover and in forward flight is further addressed in reference 6 where both 
lag-flap and flap-lag systems are studied. The arrangement of the hinges is 
shown to influence the stability boundaries for both hover and forward flight. 
In reference 7 two different sequences of rotations are used to describe the 
orientation of the cross section of the deformed elastic beam and second-degree 
nonlinear expressions are developed for torsion, the bending curvatures, and 
the [T] matrix. These expressions are used in reference 8 to derive two sets 
of equations of motion for a rotating, elastic-beam model of a helicopter blade 
in hover and in forward flight. In the two sets of equations, the torsion 
deformation is described by the appropriate angles of the two sets, which dif- 
fer from each other and from the variable of references 3-5 as well. The 
differences in the variables are not mentioned in references 7 and 8, and in 
reference 8 the fact that the formulations must be equivalent is only conceded 
to be an unproven possibility. The following conclusions are drawn in refer- 
ences 6-8: (1) the transformation sequence used in deriving the equations 

influences the stability and dynamic behavior of an elastic beam (refs. 6, 8); 
and (2) the work of references 3-5 is incorrect in some respects (refs. 7, 8). 

Because of the importance of the fundamental geometric relationships of 
slender beams to the nonlinear analysis of rotor blades, and because of the 
apparent differences that exist between references 3-5 and references 6-8, a 
detailed exposition of this subject has been undertaken in this report. The 
influence of the choice of angles is thoroughly investigated in this report. 

The fact that any of the hinged-beam approximations can be analyzed with any 
choice of angles will be demonstrated. Although the hinge arrangement is cru- 
cial to accurate analysis, the choice of angles used in modeling is immaterial 
(ref. 9) as long as singularities, discussed in section 2, are avoided 
(ref. 10). Since an elastic beam, unlike the rigid beam with its many possible 
hinge arrangements, is a single physical system, we should not expect the 
choice of angles to influence the outcome of the analysis at all. Otherwise 
a completely intrinsic formulation like the Kirchhof f— Love equations (ref. 1) 
or their recent extension by Reissner (ref. 2) would be impossible. Finally, 
the work in references 3-5, while essentially correct, does suffer from lack 
of clarity in some areas and it is hoped that this report will clarify the 
0 arlier work as well as the ensuing discussions in references 68. Part of 
the difficulty in reconciling the results of other investigations with refer- 
ences 3-5 stems from the nature of the torsion variable used therein. The 
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properties of this variable need to be clarified, especially with respect to 
its use in conjunction with Hamilton’s principle. 


1.3 Procedure 

In this report, the appropriate relations describing the geometry of the 
deformed beam are rederived and expanded to include the components of angular 
velocity and virtual rotation. The appropriate changes of variable are estab- 
lished to prove that the two formulations in references 7 and 8 are equivalent 
to each other and equivalent to the approximate form of the exact expressions 
derived herein and in references 3-5. A thorough comparison with the results 
of other work is undertaken along with a discussion of the geometric nature of 
the different torsion variables in references 3-8. The consistent reductions 
of the kinematical equations for infinite torsional rigidity and for the 
hinged, rigid-beam approximations are also examined. Finally, the application 
of Hamilton s principle to systems described by the torsion variable of refer- 
ences 3-5 is rigorously established. 


2. THE LARGE-DEFLECTION GEOMETRY OF AN EULER-BERNOULLI BEAM 


In this section the large-deflection geometry of an Euler-Bernoulli beam 
xs derived. For the purpose of discussing the geometry, we will not consider 
warp. Without warp, cross sections remain plane and normal to the deformed- 
beam elastic axis in an Euler-Bernoulli beam. Thus, the geometry is completely 
determined by three deflections u, v, w along the x, y, and z axes (fig. 2) 
and some appropriate measure of the elastic twist. We must then express all 
other geometric quantities in terms of these four deformation variables. The 
orientation of the reference axes with origin at P with respect to the 
space-fixed axes is given by elements of the [T] matrix. Also, expressions for 
the torsion and components of bending curvature, angular velocity, and vir- 
tual rotations are developed. 

We consider an initially straight beam segment with the elastic axis 
along the x axis, illustrated in figure 1. In its undeformed state the beam 
is pretwisted so that the principal axes of the cross section, denoted by y, 
and z^, are at an angle 9^(x) with the y and z axes, respectively, and 
0t(O) is denoted by 0^. The cantilever root boundary condition is considered 
here for the purpose of illustration. In principle, any set of boundary con- 
ditions could be treated by slightly modifying the method that follows. 

When the beam bends and twists to some deformed configuration, as in 
figure 2, the geometry can be completely described by the deflections u, v, 
w of the point P on the elastic axis and by the direction cosines of the' 
reference axes X3, yj, and Z3 where X3 is tangent to the elastic axis of 
the deformed beam. The matrix of direction cosines of X3, V3, Z3 with 
respect to x, y, z constitutes a coordinate transformation matrix denoted 
by [T] between the unit vectors i, j, k associated with the local reference 
axes X3, y3, Z3 of the deformed beam and the unit vectors I, J, K asso- 
ciated with the space-fixed axes x, y, z. The relationship is given by 
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( 1 ) 


i 


I 1 

j 

.= [T]. 

J J 

k 


K 1 

1 1 




When any two of the axes about which rotations occur have the possibility of 
being coincident within the primary range of interest of the rotations, singu- 
larities result. This is the disadvantage in using what are commonly known as 
classical Euler angles (ref. 10) since the first and third rotations are about 
the same axis when the second angle is zero. The remedy is to have the rota- 
tions about axes that do not approach one another for rotations in the neigh- 
borhood of zero since small deformations are a special case of interest for the 
present problem. These modified angles, as used in this development, are also 
sometimes called Euler angles in the literature (e.g., ref. 11). We assume 
that the unit vectors i, j, k are fixed in a rigid body B and are initially 
coincident with axes x, y, z. Rotations ts 6 of B occur about 

k, - j , and 1, respectively, but not necessarily in that order, so that i, j, 
k are finally aligned with the x^, y^, Z 3 axes. Therefore, the angular 
displacement 6 rotates the beam cross section about what is, for t = 3 = 0 , 
the elastic axis; the angular displacement B rotates the beam cross section 

out of what is, for t = 9 = 0 , the horizontal (x-y) plane; the angular dis- 

placement ^ rotates the beam cross section in what is, for 3=9=0, the 
horizontal plane. The first rotation always occurs about either x, y, or z. 
The orientation of the exact axis about which the second rotation occurs is 
dependent on the axis and magnitude of the first rotation. The third rotation 

always occurs about either x^, y^, or z^. In the text of this report, only 

two of the six possible sequences of angles will be considered: 3 ^, 6 ^, 

and 32 » ^ 2 * ^ 2 * subscripts are added to emphasize the fact that, for a 

given orientation of X 3 , y 3 , Z 3 with respect to x, y, z, C 2 ^ ^ 1 ’ 

and 0^ ^ 0^. Since it is clear that the choice of angles will not change the 
outcome of an analysis, it is natural to inquire at this stage why it is neces- 
sary to consider two different sequences. It is one purpose of this report to 
study the outcome of expressing the deformation kinematics in terms of the two 
sequences and to address the issues listed in section 1 . 2 . 

The axes of rotation for the 9^ sequence may be described by a 

rotation Cj about the z axis resulting in a new set of axes xj, yj, z^, 
followed by a rotation 3 ^ about the -y^ axis resulting in a new set of 
axes X 2 , y 25 finally followed by a rotation 0 ^ about the x^ axis 

resulting in the X 3 , y 3 , Z 3 axes (see fig. 3). The axes of rotation for the 

^2’ ^ 2 * ^2 sequence may be described by a rotation 33 about the -y axis 
resulting in a new set of axes x^> y^* z^, followed by a rotation ^2 about 

the z^ axis resulting in a new set of axes x^, y|» finally followed by 

rotation ©2 about the x| axis resulting in the x^, y^, Z 3 axes (see 
fig. 4). The final set of axes must be the same regardless of the sequence 
of rotations. The matrix [T] may be formulated in terms of any set of three 
angles that will rotate the axes x, y, z to X 3 , y 3 , Z 3 . The discussion 
here is limited to the sequences 9^ and 32 ? ^ 2 ^ ^ 2 * However, all 

six possible sequences of the angles 3^ and 0, as defined above, are used 

to develop the geometry in the appendix of this report. The conclusions hold 
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for any set of angles, regardless of definition, as long as they will describe 
the orientation of axes x^, y^, z^. 

The geometry of the deformed state is now developed. The position vector 
describing any point on the def ormed-beam elastic axis can be written as 
(see fig. 2) 


r= (x + u)I + vJ + 

The unit vector tangent to the elastic axis of the deformed beam is 

= (x + u)”^I + v"*”J H- w"*^ 
dr 


( 2 ) 


(3) 


where r is the curvilinear distance coordinate along the def ormed-beam elas- 
tic axis and ( = B/3r( ). Since the cross section remains normal to the 

elastic axis during deformation 


^ ” '^ 11 ^ "^ 13 ^ 


(4) 


where T^j is the element on the 1th row and in the jth column of [T]. 
Thus 


T^^ = (x + u)‘ 


Ti, = v^ 


T = w"^ 
^13 


(5) 


no matter what angles have been used to express [T]. Since the [T] matrix is 
orthonormal 


'll 


+ T 


12 


t2 

13 


= 1 


( 6 ) 


and a relationship between x and r is easily obtained from equations (5) 
and (6) 


(x + u)'^ = /I - v+2 - „+2 


(7) 


or 

r’ = /(I + u’)2 + v'2 + w’2 (8) 

where ( )' = 3/3x( ). The preceding equations must be combined with explicit 
expressions for the elements of [T] to completely specify the geometry of the 
deformed beam. 
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Explicit expressions for the [T] matrix, components of the angular veloc- 
ity and virtual rotation vectors, the torsion and components of bending curva- 
ture are now developed. The Sj, 0j sequence is chosen arbitrarily for 

this task. The [T] matrix can be determined from the successive transforma- 
tions relating x^, y^, to x, y, z using figure 3: 


[T] = 


c„ c 


'1 


?! 3^ 0^ 0J Cj 


^1 


c„ - s_ s„ s 


?1 01 


■1 ’■^1 1 


0, 6 


c„ s. 


1 '"1 




c c 

6l 01 


(9) 


T* 

where c^ ^ — cos( ), s^ ) = sin( ), and [T] ^ = [T] . The angular velocity 
is taken from figure 3: 


u) = ?iK - 3i(c J - s I) + 0ii 
9l ^1 


= w.i + (D.i + (jj, k 
1 J k 


(10) 


where ( ) - 9/3t( ). The components of u along i, j, k, determined from 
equations (1), (9), and (10), are 


""i ~ 


.. = -3iCg^ + qcg^sg^ 


‘"k = ^1^01 


i 


( 11 ) 


The virtual rotation vector Sif) can be expressed by replacement of (") with 
6( ) in equations (10) and (11) (ref. 10). Thus 


6^. = 5?iK - 33i(c J - s I) + 60 i 
^1 ^1 


6 i{;^i + + 64/j^k 


( 12 ) 


where 


6ip^ = 501 + s^ 6^1 


(13) 


10 



The virtual rotations are needed to express the virtual work of applied exter- 
nal moments (such as the aerodynamic pitching moment). The torsion (or angle 
of twist per unit length) and the bending curvatures may be deduced with the 
use of Kirchhoff's kinetic analog (ref. 1) by replacing ( ) with ( ) in 
equation (10) 


ic = - 31(0^ - s^^I) + e'J'i 


= K.i + K.j + K, k 
1 J k 


where k^, the torsion, is 


K. = 0 -J; + 


1 


(14) 


(15) 


and the components of bending curvature in the j and k directions are 




1 


’1 “1 




(16) 


We now define an angle of elastic twist (f) so that 

(0^ + (f))"^ = 


(17) 


where 0'^= 0[x+. We note that 9^ and (|) can be considered together because 
each occurs about the beam elastic axis. This angle of elastic twist is dis- 
cussed in references 12-15 and used in references 3-5 as the torsion variable. 
From equations (5) and (9) we find that 


s 


Pi 


= w+ 


(18) 


Then, we differentiate equation (5b) and make use of equations (5) , (7) , (9) , 
and (18) yielding 


+ 


+ + ++ 

V w w 


/l - V+^ - w+2 (1 - w"*"^)/! - v' 


,+2 


w 


,+2 


(19) 


Substitution of equations (18) and (19) into equations (15) and (17) yields an 
expression for in terms of variables v, w, and (|) which does not depend 

on any of the angles 6^, 0^: 


bT = (0. 


(^)+ - 


v'^w"*"w'f 


/l - v"*”^ - w'*’^ 


1 - w 


F2 


( 20 ) 
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or 


01 - + 


r 


w' 


A - v+2 


.,+2 


( 4 + , v+w+w++\ , 

+ rr^)«ir 


if OjCO) = 0t(O) = 6 and t|)(0) = 0. 


( 21 ) 


The [Tl matrix may now be expressed in terms of v, w, and with no 
dependence on Cj, gj, and ^i. We make use of equations (5), (7), and (9) 
obtaining 


[T] 


/l - v+2 _ w+2 


-v'*‘c - W+Sq /l - v+2 - w+2 

^1 ^1 

A - w"^^ 

v'’’s - w+c A - v+2 _ 

A - w+2 


c A - v+2 - w+2 _ v+W+S^ 

y 1 0 

A - w+2 

-s A - v"f ^ - w"*"^ - v+w+c 

Qj e 

A - 


S A - w'f^ 
®1 

c /I - 


( 22 ) 

where 0^ is given by equation (21). This result for [T] was first obtained 
in reference 3 in a different manner. The components of the angular velocity 
vector, the virtual rotation vector, and the rotation per^unit^length vector 
can now be expressed in terms of v, w, and 4). In doing this we make use of 
equations (5), (7), (9), (11), (13), (15), (16), and (17): 


CO . 
1 


to . 

J 


w 


/l~- v**“^ - w"^^ 


+ - V+w+w+^ ^ 


• -I- 

-w^c , 


A - w+2 


+ 


s /l-w“*"^ 
/ . 4 . v'*“w“^w'*" \ ^ 1 


• -I- 

W^S. 


A - w+2 


+ (v+ + ^ 

f2 \ 1 - w+2/ A _ 


/I - v+2 - w+2 

A - w+2 


/I — v+2 - w+2 


(23) 
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6>|^. = 60, + 

1 1 /I _ v+2 - ^2 


-C 0 5w+ 


6i);. = 

^ /l - w‘ 


w+2 \ 1 - 


/I - w+2 


/l - v”*"^ - 


s 5w+ 

6i|^, = 1 + 

^ /I - w+2 




/l - 


W / /l - - w+^ 




= (6^ + ({)) 


K . = 

^ /l - ^ 


V 


v+^ + 


v+w+w^' 


/l - w +2 


1 - w+ 2 / _ ^+2 _ ^+2 


/l - W +2 


A - W+2 


+ (v++ + 


v+w+w+^-' 


1 - vA2/ _ ^+2 _ ^2 


The variables 0, and 60, are obtained by differentiation and variation of 
0, in equation (21). Equations (21) through (25) provide a complete and 
exact description of the def ormed-beam geometry. These equations are indepen- 
dent of i;^y and 0 ^ since 0 ^ is expressed in terms of v, w, and c}> . 

To illustrate the uniqueness of the geometry we present the development 
again with a different sequence of rotations 32 > ^ 2 ’ ® 2 * make use of 

equations (1) and (5), which do not depend on the sequence of rotations. The 
[T] matrix can be determined from the successive transformation relating 
X 3 , yg, Zg to X, y, z using figure 4. 


c c 

C2 B2 


^2 ^2 


[T] = -c c s - s s 

^2 ®2 ^2 ^2 °2 




2 A ^^2 ^2 


SqC- C-^S S_ 

©2 ^2 ®2 ^2 ^2 


“Sq Cq S- S C- s_ c 

^2 ^2 ^2 ^2 ^2 ^2 ^4 




The angular velocity vector is, from figure 4 

a> = -3,J + L(c^ K - s I) + Ll 
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the components of which are 


oj. = 02 - 


"j = "^2S2 "c2 ^2^92 


(28) 


“k ■ + 62=9/5^ 


The virtual rotation components are 

6\l) = 602 - s 632 




‘♦k ■ + "o/c ,«82 ) 


(29) 


The torsion is 


Ki = e+ - 


(30) 


and the bending curvatures are 


= - 4^0 + 4^0 

j 2 ^2 2 

'k - 4 = 9 ^ + 4=9/5^ 


(31) 


Substitution of equation (17) into equation (30) yields an expression for G"*" 




(32) 


Equations (5) and (26) lead to 


s = v"*" 
^2 


4 = 


w 


,-H- 


+ 




/l - v+2 - w+2 (1 - v'''2)/i _ v+2 - „+2 


( 33 ) 


Substitution of equations (33) into equation (32) yields a relationship 
between &2 and v, w, and (() 
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^ ^ . V+W+V-H-\ 


(34) 


or 


j 


0 - = 6 . + + I - 

t I /i - v+2 - w+2 

0 


(”^ ^ r^) 


dr 


(35) 


where 62 ( 0 ) = 0^(0) = Bq and (})(0) = 0. Equations (5), (26), and (28) 
through (31) can be used to express the exact deformed-beam geometry in terms 
of V, w, and <}>. The [T] matrix is 


[T] = 


/l - v+^ - w+^ 

_V+C /l - v+^ - w+^ - w+s 

^2 ; 

/l - v+^ 

V+S /l - V+^ - W+^ - w+c^ 
^2 ^ 

/l - v+2 


V ' 


w+ 


/l - - v+w“^c . 


c /l - v+^ 

®2 


/l - v+2 

/l - v+^ - w+“ 4- v+w+s. 


-s A - 

®2 


A - v+2 


(36) 


The components of the angular velocity vector, the virtual rotation vector, 
and the rotation per~unit~length vector can now be expressed in terms of 
V, w, and (j) . In doing this we make use of equations (5), (7), (19), (26), 
and (28) through (31) : 


0 ) . = 
1 


i:* * f^) 




03 . = - W 


0 ), = 


A - v"*"^ - w3"^ 

, v+w+v+ '\ 

\ 1 - v+2/ /i _ v+2 - w+2 A - v+^ 


c A - v+^ 

®2 


S A - v+2 

®2 


(37) 


A - v+2 - w+2 A - v+2 / 
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+ y+w+6v+'^ 

i ^2 ■ ^2 V 1 - v-H^; 


&i\). = 60* - 


6i|^ . 

J 




/l - v"*"^ 6 v"^s , 


/l - v+2 - w+2 A - v+2 


( 38 ) 


= ( 


6i|». = {6w^ + 


v~^w~^6v~*' \ ^2 


/l - 6 v“^c^ 


+ 


1 v“*“ / /l - /r - 


tc . 
1 


(e^ 4^) 


+ 




c A - v"*"^ 

®2 A 

+ 


A - v3"2 - A - v+^ 


(39) 


X 


k 



yAAA^\ 

1 - v+2/ 


s A - v+^ 

A 

A - v+2 - w+2 



/l - 


/ 


where 02 and 60^ may be obtained directly from equation (35). Equations (36) 
through (39), written in terms of 02 defined in equation (35), appear to be 
different from equations (22) through (25) written in terms of 0 ^ defined in 
equation (21). It is possible to demonstrate equivalence of the two sets of 
equations by simply solving each of equations ( 21 ) and ( 35 ) for c|) and 
equating the expressions for cj). The following first integral is obtained: 


so that 


+ tan 


-1 


e 


V+W+ 


A - v"*"^ - w"*"^ 


0 


s A - v“*"^ - w"*"^ + V'^W'^C 

0^ 9 

^2 A - v“*”^ A - 

C A - V+^ - wh^ - V“*"w“*"s 

u ^ 9 

C = 

®2 A - v+2 A - W+2 


(40) 


(41) 


These relations were also obtained in reference 15. Equations (40) and (41) 
can also be derived by equating the expressions for T 23 in equations ( 22 ) 
and (36). Now, by substitution of equations (41) into equations (36) 
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through ( 39 ), the results can be shown to reduce identically to equations ( 22 ) 
through (25). Thus, the sequence of rotations used in the derivation is 
proven to be immaterial and the two sets of equations are seen to be related by 
a simple change of variable. The fact is further illustrated in the appendix, 
in which all six possible sequences are examined and shown to be equivalent. 


3. COMPARISON WITH OTHER WORK 


The results of the development in section 2, except for slight differences 
in notation, are identical to the results in references 3-5. Several discus- 
sions of the material in references 3-5 appear in the literature and it is the 
purpose of this section to examine these discussions. 

In references 6 - 8 , the [T] matrix, the torsion, and the components of 
bending curvature are obtained for both the Ci> 3^, sequence and the 

82 , ? 2 ’ sequence of rotations. The expressions therein are approximated 

to second-degree accuracy in the deflections u, v, w and or 02 * The 

authors denote 0 ;^ and 0 2 by the same symbol even though 0 j and 02 are 
different angles. After recognizing this slight inconsistency in notation, 
however, it is quite easy to verify that the geometric quantities expressed 
tjierein are equivalent to the two formulations in the preceding section. 

Indeed, the development in the preceding section shows that the sequence of 
rotations used in describing the orientation of the local principal axes of 
the beam cross section is arbitrary. Thus the sequence of transformations 
cannot influence the outcome of a stability or response analysis. The equa- 
tions appear to be different but a change of variable can always be found that 
will transform the equations derived based on one sequence to that of any 
other. The application of this fundamental idea to the rigid-beam approxima- 
tion of an elastic beam will be discussed in section 6 . 

Let us now examine the question of whether the sequence of transformations 
used in an elastic-beam analysis can affect the results of the analysis as 
discussed in reference 6 . This involves demonstrating that the two develop- 
ments in references 6-8 are actually equivalent to each other and to that of 
references 3-5 and that the [T] matrix and kinematic quantities developed in 
references 3 — 5 , despite the apparent differences with the results of refer 
ences 6 - 8 , are nonetheless correct. It is convenient for comparison with 
references 6-8 to approximate the geometric quantities obtained in section 2 
to second order in bending and torsion deflections. 

Let us consider the case with 0^ = 0 and v+, vT*", 0^, 62 and of 
0(e) so that e^ << 1. Equation ( 8 ) shows that 


,i2 = 


= (1 + u')- 


+ v' ^ + w 


'2 = 1 + 2 e 


(42) 


where e is the component of longitudinal strain of the elastic axis from 
Green's strain tensor. In engineering beam theory, e is neglected with 
respect to unity. Here we assume e = O(e^) so that r' = 1 + 0(e ) and 
u' = O(e^); thus, derivatives with respect to r may be replaced by 
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derivatives with respect to x. For comparison we write only the Tj, and T,, 
elements of the [T] matrix, the torsion, and the bending curvatures of the 
? 1 » formulation: 


^23 “ 

'^32 " + 0 (e 3 ) 

= 0j + v"w' + O(e^) ► 

Kj = -w" + e^v" + 0 (e3) 


(43) 


where equations (21), (22), and (25) have been used. For the , c , 6 
formulation these quantities are * ^ 


T23 ~ ^2 ~ v'w' + O(e^) 

T32 = “02 O(e^) 

<1 = 02 - v'w" + 0 (e 3 ) > 

= -w" + e^v" + 0(e3) 

= v" + 02W" + 0(e3) 


(44) 


where equations (35), (36), and (39) have been used. These quantities are 
identical to the corresponding expressions obtained in references 7 and 8 . 
Equation (40) yields, to second order, the following relationship between 
0 j^ and e^: 


02 = 0^ + v'w' + O(e^) ( 45 ) 

which, upon substitution into equations (44), produces equations (43) identi- 
cally. These expressions, which form the basis for deriving the equations of 
motion, are thus equivalent. Equation (45) can be thought of as a change of 
variable. If it is applied to the equations of motion based on the 4 9 

sequence (flap-lag-pitch) derived in reference 8 , the result will be identi- ^ 
^^lly the other set of equations of motion in reference 8 that are based on 
the 4 ^, 3 ^, 6 ^ sequence (lag-flap-pitch). 

Equations (43) and (44) both reduce to the same expressions when the 
elastic angle of twist is introduced, from equations ( 20 ) and ( 21 ) 
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( 46 ) 


(f)' = + v"w' + O(e^) = ^2 “ O(e^) 

J *x 

v'Vdx + O(e^) = ^2 “ J ■•■ O(e^) 

0 *'o 


or 


+ 



v"w’dx + 0 (e3) 


v'w"dx + O(e^) 


(47) 


Now equations (43) and (44) become 

/ X 

v"w'dx + O(e^) 

0 


• X 

“ J v'w"dx + O(e^) 

•'o 

, I C48) 

= <t>' + O(e^) 

K . = -w" + (t)v" + O(e^) 

: 

= v" + (|)w" + O(e^) j 

Equations (47) relate <t>, 0,, and 0 2 and may also be regarded as changes of 

variable relating the angles Oj and 02 to the elastic twist angle When 

these changes of variable are applied to the development in reference 8, the 
two sets of equations derived therein based on Cj. and based on 

Bo, Co, Bo (when specialized to hovering flight) are found to be equivalent to 
the set of equations found in references 3-5 (which considered only the hover- 
ing flight case). Because of the nature of (|), this change of variable is 
most easily applied to the energy and virtual work expressions, however, not 
directly to the partial differential equations, as shown in section 5. The 
outcome of the analysis does not depend on the sequence of rotations and the 
work of references 3—5 stands as correct. 

other authors (refs. 16 and 17) have had difficulty reconciling their 
work with references 3-5. In reference 16 it is stated that the results of 
references 3-5 are in error because of an incorrect expression for torsion, a 
comment similar to the conclusions of references 6-8. The [T] matrix in 
references 3-5 is stated to be incorrect in references 16 and 17 because the 
pretwist angle is grouped with the torsion variable ratber than being imposed 
before the elastic deformation. Furthermore, it is stated in reference 16 
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that the quantity 9,- must appear in [T] . The pretwist 0. can be grouped 
with 0 ^, © 2 , and 0 in the present formulation because all of these rotations 
occur about the same axis. The development in section 2 of this report is 
exact for an Euler-Bernoulli beam, and 6 [ terms do not (and cannot, because 
of dimensional reasons) appear in [T]. In contrast to references 16 and 17 
there is no difference in the final form of the [T] matrix that results from 
interchanging the order of pretwist and elastic deformations (see eq. (17)) 

In both references 16 and 17, the two sets of axes that the [T] matrix relates 
are not the same two sets of axes that the [T] matrix in this report (and in 
refs. 3-5) relates. In references 16 and 17, the [T] matrix relates the local 
beam cross-section principal axes after deformation to the local beam cross- 
section principal axes before deformation and not to the x, y, z axes used 
here. Thus, for 6 ^. 0 a direct comparison of the [T] matrix and other ' 

geometric quantities in references 16 and 17 with those of references 3-5 and 
the present report is not possible. 

In reference 12, two of the present authors established that the lateral 
buckling load for uniform, slender cantilever beams was identical for two sets 
of equations — one in terms of 9 ^ and the other in terms of ip. These 
equations were different in appearance, yet they were equivalent and yielded 
the Identical buckling load. The development in reference 12 in terms of 
proceeds from the Klrchhof f-Love equations (ref. 1 ) which in no way depend on 
the sequence of rotations used to describe the cross-section orientation. 
Careful examination of the different approaches noted in reference 12 illus- 
trates that it is important to understand that the two geometric angles, 
denoted by^ 0 ^ and 6 ^ in this report, and the elastic angle of twist ({> are, 
in fact, different quantities. Some of the confusion about the influence of 
different sequences of rotations evident in references 6-8 apparently stems 

from their use of the same symbol for the three different quantities 0, 0 

and (() . 1 ’ 2 ’ 


In reference 18, nonlinear equations of motion for an Euler-Bernoulli 
beam are derived for the purpose of investigating stability of the beam with 
applied loading. Classical Euler angles are used to describe the cross- 
section orientation. The three axes about which rotations occur for infini- 
tesimal angles are not mutually orthogonal for classical Euler angles, for 
which two rotations occur about the axis that lies along what is nominally the 
beam elastic axis. Thus, the torsion variable x is defined in reference 18 
as the sum of the first and third Euler angles. The torsion kj and the I 
component of the virtual rotation vector 6 ijj = 6* • I are given to second 
order in reference 18 as ^ 




T ' + Y (v"w' - v'w") + 0 (c 3 ) 

6 t - Y (<5v’w' - v' 6 w') + 0 (g3) j 


(A9) 


9 • • 

This is not true of the other four sequences treated in the appendix, 
however. There, angular displacements 63 , 0 (^, and 0 . do not rotate the 
beam cross section about the deformed-beam elastic axis. 
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These quantities, in terms of 0j^ and 02> are, from equations (12), (25), 
and (40) 


K. = e! + v"w' + O(e^) = - v'w" + 0(6^) 

1 1 

64)^ = 60J + v'5w' + O(e^) = 602 - w'6v' + 0 ( 0 ^) 


( 50 ) 


Equations (49) and (50) can be shown to be equivalent to second order by use 
of the following changes of variable: 


T = 0 , 


+ 0(e3) = 0„ - 


O(e^) 


(51) 


Thus, the equations of motion derived with the aid of classical Euler angles 
in reference 18 can be shown to be equivalent to those derived with the rota- 
tions defined in this report. To formulate the complete, exact geometry using 
classical Euler angles and the elastic angle of twist, is quite difficult, 
however, and may not even be possible because singularities appear in the 
integrals when deformations are small instead of large as in the present for- 
mulation. Also, the angles used in this report, unlike classical Euler angles, 
bear a strong resemblance to the geometric change in orientation during 
moderate deformations. 

The results of this section show that the expressions for def ormed-beam 
geometry developed in section 2 and in references 3-5 are equivalent to a 
variety of formulations found in the literature. The results of references 3-5 
are also confirmed. The nature of the torsion variable ()) and how a formula 
tion in terms of 4 differs from one in terms of 0^ or Op has not been 
treated in the literature and will be addressed below. 


4. NATURE OF THE TORSION VARIABLES 0^, Op, AND 4> 


In this section we will examine the nature of the torsion variables 
0,, 0 , and (t>. The kinematical relationships that define these variables are 
used to establish differences among them. An example problem is Introduced to 
aid in illustrating these differences. The geometric relationship between 
0, and 0p is examined by projecting the local cross-section principal axes on 
the yz plane. Finally, the elastic twist 4 > is shown to be a quasi- 
coordinate (ref. 10) in the space domain because it is related to physical 
angles through Integrals, and its mathematical analogy with the axial deflec- 
tion of a beam due to longitudinal strain is established. 


4.1 Kinematical Differences 

The elastic angle of twist 4> is defined to be the elastic component of 
the integral of the torsion Thus, from equation (17) 
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( 52 ) 


<!> 



K. dr - 0 ^ 
1 t 


The elastic torsion moment is then simply 


GJ(j)‘^ (ref. 4). 


The angle may be thought of as the beam twisting caused by a distrib- 
uted, externally applied torsion moment acting about the local x- axis. The 
angle (ji does not completely define the orientation of the reference axes at 
point P even when and 6 ^ (or v+ and w+) are known at P. In addition, 6 , 

(or 62 ) must be calculated to obtain the complete description of the geometry 
of a deformed beam. This requires knowledge of and for all x inboard 

of P. Both 0^ and ©2 occur as rotations about X 3 , the deformed-beam elas- 
tic axis (i.e., the i unit vector). By definition, 0^ and ©2 are the third 
angles in the two different sequences gj, 0 ^, and g 2 , ? 2 » ^ 2 » respec- 

tively, as described in section 2. Since the final orientation is fixed, the 
angles 0^ and ©2 inust differ. The fact that the elastic twist is not, in 
general, equal to either of the angles 0 ^ and ©2 is simply because of the 
kinematic effects of combined bending of the beam in two directions, v and w. 

If the bending deflections and elastic twist angle are known, 0 , and 6 , can 
be written according to equations ( 21 ) and ( 35 ). 


0 , = 0 


J r 

L-h- + 

/— - v+2 _ w+2 V 1 - w+V 


dr 




J /l - v +2 - w +2 \ 1 - V+ 2 / 


(53) 


They can be approximated quite well by the simpler expressions: 


'1 = 6^ + <)> - r v"w' 


dx + 0 (e 


^2 “ 


+ d) + 


r v*w"dx -h O(e^) 


(54) 


Thus, rotations about the beam elastic axis (geometrically like elastic tor- 
sion rotations) may be produced by combined bending deflections alone, even 
without elastic twist (p or pretwist 0^. 


4.2 Example Problem 

The following example problem is included to further confirm and illus- 
trate the differences between 0 ^ and 0 ^ as well as to demonstrate the use of 
the geometric relations developed in section 2. Consider a beam that is first 
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inclined by an angle 3 q from the x-y plane and then bent into a circular 
arc, without torsion or elongation, in the plane of its elastic axis and major 
principal axis as shown in figures 5 and 6. Thus 


= 0 


Kj = 0 


<k = -Kq = constant J 

The equations for the space curve of the beam elastic axis are, from the geom- 
etry of figures 5 and 6, 


cos sinCK^r) 


1 - cos(k r) 


- r ; V = 


sin 3 q sinCK^r) 





Figure 5.- Deformed beam for example problem. 



Figure 6.- View of deformed-beam elastic axis in the plane determined by 

i and j (view AA from fig. 5) . 


The matrix of direction cosines at any location is quite easily obtained from 
figures 5 and 6 



We wish to show that given u, v, w and for a beam, the equations of 

section 2 will provide an exact description for the deformed-beam geometry. 
Since is zero then (j) = 0. The geometric angle 6 ^ can be calculated 

directly from substitution of equations (56) into equation ( 21 ) yielding 

= tan"^ [sin(KQr)tan e^] (58) 


Substitution of equations (58) and (56) into equation (22) for [T] yields 


[T] = 


'cos(K:Qr) cos 3 q -sin(KQr) cos(KQr)sin gj" 


sin(K .r) cos 


'0^ 

-sin 


cos(iCQr) sin(Kp|r)sin g 
0 


^0 

cos g 


0 


(59) 


which is identical to the matrix of direction cosines defining the example 
deformed beam in equation (57) . It can be verified by substitution of equa- 
tions (56) and (58) into the bending curvature formulas, equations (25), that 
the exact relations defining the problem in equations (55) are recovered. 

Thus the formulas of section 2 provide means of exactly determining the orien- 
tation of the deformed-beam cross section when the elastic bending and elastic 
torsion deflections are specified. 


The same results for the transformation matrix and bending curvatures are 
obtained if one uses the expressions derived from the g 2 , notation 

sequence, equations (35), (36), and (39). Proceeding as above but starting 
with equation (35) we find 


02 = 0 


(60) 


and from equation (36) 


[T] = 


cos(KQr)cos 


-sin(Kpr) cos(K^r)sin 


sin(K r)cos g^ cos(Kp,r) sin(K.r)sin g. 


-sin 


(61) 


cos g 


0 


The transformation matrix is the same as obtained in equation (59) and thus 
again duplicates the example problem transformation matrix given in equa- 
tion (57). The fact that 9*^ / 02 should clearly indicate that the differ- 
ences between 9^ and 62 cannot be ignored. If 0^ and 0^ were treated as 
the same angle it would not be possible to correctly describe the geometry of 
the example problem in terms of both 0 ^ and 02 • 


4.3 Geometric Relationship Between 0^ and 02 

It is helpful to visualize 0^ and 02 for the general case of coupled 
bending deflections. This may be accomplished by considering the local 
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Figure 7.- Projection of the local 

cross-section principal axes on the 
y-z plane. 



Figure 8 .- The projection of unit vec- 
tors 2 ^ on the y-z plane. 


cross-section principal axes pro- 
jected on the y-z plane. In fig- 
ure 7 the projection of the y^ and 
z^ axes (the unit vectors j and k) 
are shown in the y-z plane (the 
plane containing unit vectors J and 
K). The beam cross section appears 
skewed because it is not in the plane 
of the J and K unit vectors. The 
projections of j and k, denoted by 
jp and kp, are given by 






(62) 


By inspection from figure 8 


sin 


sin 




'23 


+ Ti 


V 22 23 

-T. 


/ 


32 

T’2 ^ rp2 

^32 ^ ^33 


(63) 


Expressions for the matrix components 
T 23 and T ^2 terms of 0 ^ or 62 

can be found in equations ( 22 ) 
or (36) , respectively 


T ,3 - s A - 



(64) 


Combining equations (63) and (64) and dropping higher-order terms, we find 
that 


= 9^ + O(e^) 

02 = 02 + 0(e3) 

and, according to equation (45) 

02 - Of = v'w' + 0 (e 3 ) 


(65) 


( 66 ) 
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These results show that for small angles, Qi and ©2 are essentially the kine- 
matic rotations of the j and k axes, respectively, about the i axis. The 
fact that they are not, in general, equal is because the 3 -k plane does not 
remain parallel to the J-K plane. 

Returning to figure 5, the view of the beam in the direction of the I 
axis provides a clear illustration of the kinematic rotation of the beam. The 
orientation of the projections of the j and k axes in figure 5 is consistent 
with the geometric interpretation of 0^ and 02 in equations (65) and (66) and 
with the example calculations showing that 0^ = -v’w’ and 02 == 0. These 
kinematic rotations, devoid of elastic twist, were previously identified in 
reference 3 where their effect on the blade pitch angle was noted. Hence, the 
term "kinematic pitch rotation" was used to describe these effects. 


4.4 Mathematical Analogy of the Quasi-Coordinates and u^ 


It is well known that even in the absence of longitudinal strain of the 
elastic axis, there may be significant nonzero axial deflection of a beam in 
bending. The above calculations show that even in the absence of torsional 
(shear) strain, there may be some nonzero pitch rotation of the cross section 
of a bent beam. To examine this analogy in detail we must first develop the 
appropriate relations concerning axial deflection. The longitudinal strain of 
the elastic axis may be explicitly determined from the relationship of x and 
r given in equation (42) . Thus 


e 


= u’ + 





(67) 


We now let u = u^ + u^ where is a kinematic axial deflection due to 

bending only and u^ is the axial deflection due to longitudinal strain of 
the elastic axis. The strain becomes 


e . u; + u; + 


<u; + 


r» 2 


2 


The strain must vanish when u^ = 0. Thus 


,»2 


% 


+ 


r'2 


+ 


w 


I 2 


= 0 


( 68 ) 


(69) 


or 

u^ = A - - w’^ - 1 (70) 

The kinematic axial deflection, sometimes called foreshortening (refs. 19-21), 
is then given by 
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( 71 ) 


= J /l - dx - X 

•'o 

= - j (v'2 + w'2) + 0(e‘*) 


Thus 


/ X 

/l - v’ ^ - w’ ^ dx 

0 

= u + y r + w*^)dx + O(e^) (72) 

•^0 

Comparison of equations (72) and (46) reveals that u^ and <|) are mathematical 
analogs. The deflection u^ is devoid of kinematic shortening due to bending. 
The elastic twist (J) is devoid of kinematic pitch rotation due to bending. In 
the time domain, quasi-coordinates (ref. 10) are usually based on integrals of 
either components of velocity or angular velocity. By Kirchhoff’s kinetic 
analog (ref. 1) the angular velocity components carry over to bending curva- 
ture and torsion relations while the velocity component along the elastic axis 
can be shown to correspond to a longitudinal strain approximately equal to e. 
To obtain the latter correspondence, the velocity of point P must be written 
subject to the constraint that the velocity vector is tangent to the beam 
elastic axis. The component of strain from the analog is the difference of 
the velocities, for (*) = ( )"^, with and without elongation of the elastic 
axis. Hence, both u^ and cj) are denoted as quasi-coordinates in the space 
domain because they are related to the physical displacement u and angles 
0^ and 02 through integrals and because their derivatives have velocity and 
angular velocity components as kinetic analogs. 


5. QUASI-COORDINATES IN THE SPACE DOMAIN AND 
HAMILTON’S PRINCIPLE 

The angle of elastic twist (j) has been shown in section 4 to be a quasi- 
coordinate in the space domain similar to the elastic component of axial 
deflection u^. Although several derivations appear in the literature, making 
use of at least one of these variables (refs. 4-6, 8 and 19-21), it is well 
known that the standard form of Lagrange’s equations (and thus Hamilton’s 
principle) is not adequate when quasi-coordinates in the time domain are used 
(e.g., see ref. 10). It is natural to ask if there are similar obstacles when 
the above quasi-coordinates are used in the standard form of Hamilton’s 
principle. To answer this question we will consider both a general development 
and a simple example derivation. 
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5.1 General Development 


We will first write symbolically the standard extended form of Hamilton’s 
principle for nonconservative systems using u, v, w, and For simplicity, 

we will consider only a linearized stability analysis so that time dependence 
is factored out of every term. We then will rigorously take the variation. 

The change of variable will be made to variables u^, v, w, and after the 
variation. The resulting expression, if there are no obstacles, should corre- 
spond identically to the expression obtained by starting the formulation in 
terms of u^, v, w, and c|) and then taking the variation. 

The standard extended form of Hamilton’s principle for an unpretwisted , 
nonwarping Euler-Bernoulli beam with variables u, v, w, and 9^ is given by 



where H is the kinetic energy minus the potential energy and (W is the 
virtual work of the nonconservative aerodynamic forces. To simplify the fol- 
lowing development, equation (73) is written for a linearized stability analy- 
sis, with time dependence factored out of every term or a static analysis 
without time dependence. The variation of H when expanded yields 
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where the virtual work has been expressed in terms of the generalized applied 
loads and = 60^ +w’6v’ to second order. The forces are resolved along 

the x, y, z axes, respectively, and the pitching moment acts about the x^ 
axis. We now make the change of variable that transforms equation (74) to a 
different expression which involves u^ and c(). The relationships between 
u, 9^, and u^, 4) are, to second order, from equations (46) and (72) 
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Differentiation of equations (75) yields 
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It should be noted 
have been made if 
as might occur for 
into equation ( 74 ) 


that the above variation in terms of (or u^) could not 
(or U 3 ) were also an explicit function of 9 ^ (or u) , 
general quasi-coordinates. Substitution of equation (76) 
yields 
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+ (60^ + w’ 6 v') dx = 0 (78) 




Up to now the change of variables in Hamilton’s principle is equivalent to a 
direct substitution for u^ and cj) into the partial differential equations 
corresponding to equation (74) . The single underlined terms appear in the 
V and w equations from u^, and the double underlined terms are from Not 

all of these terms are necessarily visible at this step, however, depending on 
the order of terms that are retained. The example given below should clarify 
the nature of these terms. We now substitute for 60^, 69|, 6 u, and 6 u’ in 
equation (78) 
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All of the underlined terms cancel out leaving 
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This is precisely the same expression that would have resulted from starting 
with the quasi-coordinates and taking the variation directly. Since the vir- 
tual work term involves only a straightforward change of variable, Hamilton's 
principle is valid when written in terms of quasi-coordinates Ug and 4 ), Thus 
the application of Hamilton's principle with variables Ug and 4 reduces to 
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where the following identity was used: 

J* f(x) r g(xj)dx^ dx = r g(x) r f(x^)dx^ dx 

•'0 ^0 ^0 


(82) 


It is observed that the use of the quasi-coordinates Ug and 4 has the effect 
of replacing the terms underlined in equation (78) with integrals of the 
applied loads from both the u and 0 ^ equations. 


5.2 Example Derivation 

An actual derivation of equations will be helpful in understanding the 
results from the general development above. Consider a rotating cantilever 
beam that undergoes coupled axial deflection and flap (out-of-plane) bending 
only. For simplicity we consider only the static equilibrium equations (no 
time dependence) for a beam with a vertical distributed load proportional to 
X + u and an axial distributed load P. The beam is allowed to be extensional 
in that elongation of the elastic axis is allowed. The important special 
case where the elastic axis is inextensional is also considered. The position 
vector for a point at a distance C from the elastic axis is 


r = (x + u)I + wK + 5 k (83) 

and the axial strain, neglecting elongations with respect to unity (ref. 22 ), 
is 
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Application of Hamilton’s principle yields 



where K 
equations 


+ + w'6w') + EIw"[(l + w'2)6w" + w"w'6w’] 

- P6u - m^2^(x + u)6u - K(x + u)6w|dx = o (85) 

is a given constant, A = J*J dA and I = // C^dA- The Euler-Lagrange 
and boundary conditions follow immediately: 
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Let us consider for the moment the case of an inextensional beam. To special- 
ize the equations we must do two things: (1) solve equation (86a) for 

EA[u’ + (w’^/2)] and substitute the result in equation (86b); and then 
(2) solve for u in terms of w and substitute into equation (86b). Although 
u’ + (w’^/2) = 0 for inextensibility, EA ^ for this case and the product 
EA[u* + (w’^/2)] is finite, and nonzero in general. For step (1) 
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Thus, the single w equation for an inextensional beam becomes 



with boundary conditions as above. 

We now return to the extensional case and make use of the quasi-coordinate 
Ug so that 
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with boundary conditions as_above. The underlined term in equations (92) 
and (93b) corresponds to (9H/3ug)w' 6w* in the Hamiltonian (see equation (78)). 
Equations (93a) and (93b) are identical to the result of substituting equa- 
tions (91) into equations (86a) and (86b). For the case of inextensibility, 
we still must follow two steps: (1) solve equation (93a) for EAu’ and sub- 

stitute into equation (93b) ; and (2) substitute Ug = 0 into equation (93b) . 
The result is equation (90) . 
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To express Hamilton’s principle in terms of and w alone requires 
substitution of 
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The Euler-Lagrange equations are 
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For the case of inextensibility, we simply set u^ =0 in equation (96b) and 
obtain equation (90) directly. 


Although no EA terms appear in equation (96b), it is still suitable for 
use, along with equation (96a), when EA is finite. The reason for this is 
now clear: when the analysis is formulated in terms of the quasi-coordinate 

Ug at the outset, the substitution for the tension EAug is automatically 
taken care of. The rest of equation (96a) is, in effect, integrated and sub- 
stituted automatically. This example explains why no EAu^ terms appear in 
the V and w equations of references 19-21. They do appear in the v and w 
equations of reference 8 because the case of the tension axis being offset from 
the elastic axis is considered. 


At this point, it seems appropriate to comment on the advantages of using 
the quasi-coordinates u^ and c() . Both variables tend in some cases to simplify 
the derivations of actual equations of motion (refs. 3-5, 8, and 19-21). It 
is also evident that taking the limit for infinite axial or torsional rigidity 
may be somewhat simpler. 
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6. IN THE LIMIT OF INFINITE TORSIONAL RIGIDITY 


Another way of looking at the question of the equivalence of the equations 
in references 3-5 and the two sets of equations in reference 8, demonstrated in 
section 3, is to consider the special limiting case of infinite torsional 
rigidity. It is clear that no matter what variables are used in the analysis 
for torsion there is a unique set of equations governing bending deflections 
V and w when GJ tends to infinity. For that case, all GJ terms must 
vanish from the v and w equations in order to preserve the proper mathemati- 
cal structure of the equations. This is similar to the result of the last 
section in which all EA terms vanished for the special case of infinite EA. 

For equations written in terms of ((), such as those of references 3-5, 
there are no GJ<^' terms in the v and w equations just as there are no 
EAUg terms in the v and w equations of references 19-21. With the order- 
ing scheme of references 3-5, all of the terms from the (j) equation that are 
integrated and substituted into the v and w equations automatically, through 
the process of using the quasi-coordinate (f), turn out to be negligible (i.e., 
as stated in reference 3, the torsion moment is one order of e smaller than 
the bending moments). For infinite GJ, the angle of twist per unit length is 
simply the pretwist component 0^; the elastic twist (f) and its variation 6(j) 
are simply set equal to zero because elastic twist about the elastic axis 
cannot occur. 


For equations written in terms of 0^ or the infinite torsional 
rigidity constraint is more involved. If 0^ or 02 were simply set to zero 
nonzero GJ coefficients would remain and the equations would break down. 
Instead, the GJ coefficient, equal to - 0^, must be set equal to zero 

and used to eliminate 0^ or 02* The following kinematic relations hold to 
second order (see eqs. (46)): 
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Then, the elastic torsion moment GJ(k^^ - 0^-) which is finite and not neces- 
sarily zero, must be determined from the torsion equation and substituted into 
the V and w equations. In this way all 0^ or ©2 terms and all GJ terms 
are eliminated from the v and w equations. Finally, if the aerodynamic 
pitching moment is not neglected, the virtual rotation in terms of v and w 
alone must be determined from equations (24) or (38) and (97): 
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for either 6j or The total effect of pitching moment is then in the terms 

-(w"f^Mx dxj)' in the v equation and dXj)' in the w equation. 

3 X 3 

It may be shown that these operations will yield a single set of equations in 
V and w regardless of which /ariable 9 ^ or 02 is used. The only second- 
order term is the integral term correction to the pitch angle, the "kinematic 
pitch rotation" term, first identified in reference 3. The other places that 
the integrals appear are in third— and higher-order terms and not particularly 
significant. Therefore, it is not correct to simply set 0^ or 02 equal to 
zero for general nonlinear analysis of torsionally rigid beams. In fact, when 
any physical displacement variable u, v, w, 0^, 02 that is part of nonlinear 
Euler-Bernoulli beam equations must be eliminated because of infinite axial, 
bending or torsion stiffness, it must be done by a process similar to the one 
outlined above. We have shown that u 0 when EA and 0^ ^ ©2 ^ 0 

when GJ It is just as easily shown that neither v nor w vanishes when 

the bending stiffness in that direction tends to infinity. A single set of 
equations, consistent to second order, for a torsionally rigid rotating beam 
can be obtained by simply setting <j) = 0 in references 4 and 5 or by going 
through the process outlined above for either set of the equations in refer- 
ence 8. 

Another result of equations (97) is that three angles are required to 
describe the general orientation of a reference frame in space. Thus, even 
for infinite torsional rigidity it is possible for beam bending deformations 
to place beam elements in such a position that transformations based on only 
two angles cannot describe the orientation of a beam element. This is seen in 
the example problem of section 4, in which three nonzero angles may be required 
to describe the orientation of a beam element, despite the fact that = 0. 

The required third angles are given to second order by equations (97) . 


7. RIGID BEAM APPROXIMATIONS 


When cantilevered elastic beams are modeled by an approximate system 
that substitutes a rigid, hinged, spring-restrained rod for the elastic beam, 
three degrees of freedom are often considered: an angular deflection in 

the x-y plane (when 3 = 0 ” 0) or lead-lag; 3, an angular deflection out 
of the x-y plane (when ; = 0 = 0) or flap; and 0, the pitch angle (when 
^ 3 = 0). The nonlinear equations for this case are especially interesting 

because six different hinge arrangements can be defined which will lead to six 
different sets of nonlinear equations (the linear equations being the same) , 
Simpler models, assuming no pitch degree of freedom, have been used in pre- 
liminary investigations of helicopter rotor blade stability (e.g., ref. 23). 
For this case there are two possibilities: flap-lag, where the lead-lag hinge 

flaps with the blade, and lag-flap, where the flap hinge lead-lags with the 
blade. The two hinge arrangements define two different physical systems with 
different dynamic properties although the differences are not large. Differ- 
ences in the stability of these two systems were identified and discussed in 
reference 3, and reference 6 added further numerical results and extensive 
discussion. It was concluded in reference 6 that two different sets of 
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Equations for a single Bl'CLStic beam would result from different rotation 
sequences and that this could be expected to influence the stability of the 
elastic-beam model. This incorrect expectation was based on the correct obser- 
vation that rigid-beam models with two different hinge arrangements exhibited 
differences in stability. In fact, however, a single set of equations for an 
elastic beam can be used to develop the exact equations for a rigid-beam model 
with any hinge arrangement. In so doing one must carefully consider the 
dependence of the v and w displacements on the hinge arrangement. Further- 
more it can be shown that the sequence of rotations used to develop the equa- 
tions is arbitrary for any rigid-beam hinge arrangement as well. For simplic- 
ity we will perform these operations for a second-order analysis. 

Let us assume first a lag-flap hinge arrangement with the lead-lag hinge 
at the root and the flap hinge offset from the root at axial distance e. 

The limit must be taken as c ^ 0 to recover the correct hinge arrangement 
even for the coincident hinge model used in reference 6. The following dis- 
placement functions are assumed 


v=Cx, v’ = i;, v" = 0 x>0 
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arrangement, we need to reverse the expressions in equations (99) so that 
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To reduce the elastic-beam equations in references 4 and 5 to the torsionally 
rigid case we simply set 4> = 0. The resulting equations would then contain 
the integrals in equations (97) that have different values for the two hinge 
arrangements defined by equations (99) and (100) . Evaluation of the integrals 
for the lag-flap and flap-lag hinge arrangements yields 
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for lag-flap: 


for flap-lag: 



( 101 ) 


( 102 ) 


The differences between equations (101) and (102) result in two different sets 
of equations for the two different hinge arrangements, and these equations 
are, to second order, identical to the rigid-blade equations derived in 
reference 6 . 


For two degrees of freedom 3 a-^d c, and a given hinge arrangement, any 
sequence of rotations (e.g., Ci» 3i, 9^ or 62 , ^ 2 ^ ^2^ formu- 

late the equations although it would certainly be more natural to use the 
angles associated with rotations about the hinges. The reconciliation of 
three angles with only two degrees of freedom is made by use of a holonomic 
constraint equation that relates 0 to 3 and C (ref. 10, p. 55). For 
example, consider writing equations for a flap-lag hinge arrangement in terms 
of 6 ^, Oj, instead of in terms 32 ^ ^2 ^^2 needed because the 

axes about which 02 and ^2 occur are along the hinge axes) . Thus 02=0 
and equation (40) yields 



or 


0 ^ = -tan ^(sin 3^ ~ ^ O(c^) (104) 

In fact, as long as the proper constraint is accounted for, any set of vari- 
ables can be used to write the equations for any physical hinge arrangement. 
When the constraints are properly accounted for, the equations of motion for a 
given physical hinge arrangement are equivalent, although different in appear- 
ance, for any sequence of rotations used to describe the beam orientation. 

The dynamic behavior of the system is then independent of the transformation 
sequence used in describing the orientation of the blade. The equations and 
dynamic behavior however, depend on the physical arrangement of the hinges. 
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8. CONCLUSIONS 


In this report the large deformation geometry for an Euler-BernouIIi beam 
has been developed. Some aspects of the beam kinematics discussed in refer- 
ences 3-8, 12, 16, and 17 have been clarified. The following points summarize 
the results of this study: 

1. The large-deformation kinematics for an Euler-Bernoulli beam are 
developed, including the transformation matrix relating the local principal 
axes in the deformed state to space-fixed Cartesian axes, the components of 
angular velocity and virtual rotation vectors, the torsion, and the components 
of bending curvature. The values of all the geometric quantities are unique 
at a given instant in time, but the form of the expressions themselves may 
depend on the variables used in describing the orientation of the cross sec- 
tion during deformation. 

2. The angles 3^, 6 ^ are not the same as the angles 62 » C 21 ^ 2 ’ 

both sets of which are discussed in the text. Although the geometry is unique, 
the use of the different sequences of rotations may change the appearance of 
the geometric expressions. 

3. The exact expressions for the large-deformation geometry that are 
derived in this report do not depend on the sequence of rotations used in 
defining the beam cross-section orientation. Some of these quantities were 
also derived previously in references 3 - 5 . 

4. The stability and dynamic behavior of an elastic beam are not depen- 
dent on the choice of angles used to describe the orientation of the local 
beam cross section during deformation. In fact, the two sets of equations 
derived in reference 8 , based on two different sets of angles, are shown to be 
equivalent to each other and to the equations of references 4 and 5 . 

5. The exact relationship of the pitch rotation angles 0^ and 0^ and 
the elastic twist angle (|) , established in reference 3 , is verified along with 
the exact transformation matrix first derived therein. The integrals asso- 
ciated with these quantities are correct as originally derived in reference 3 , 
and serve to show that there can be a kinematic pitch rotation of the beam 
cross section that depends on bending alone, even when there is no pretwist 
and no elastic twist. 

6 . The kinematic pitch rotation due to bending only that may occur even 

in torsionally rigid beams is seen to be mathematically analogous to a kine- 
matic shortening that may occur even in axially rigid beams due to bending 
only. The torsion variable used in this report is the elastic angle of 

twist, which is devoid of kinematic pitch rotation due to bending alone. Its 
mathematical analog is the axial deflection due to longitudinal strain of the 
elastic axis u^, which is devoid of kinematic axial deflection (foreshorten- 
ing) due to bending only. 

7. Both the elastic twist angle and the axial deflection due to 
longitudinal strain of the elastic axis u^, because they are related to 
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physical deflections through integrals and have integrals of angular velocity 
and velocity components for kinetic analogs, belong to a certain class of 
variables called quasi-coordinates. They have both been used in Hamilton's 
principle without the rigorous justification provided in this report. Although 
not previously established, their use in references 4, 5, 8, and 19-21 is 
found to be valid. 

8. Three angles are shown to be required in general to describe the local 
cross-section orientation for an elastic beam, even if the torsion stiffness 

is infinite. Only one set of equations exist for bending deflections v and 
w of an inextensional , torsionally' rigid Euler-Bernoulli beam, regardless of 
the three angles used in the derivation. 

9. When the proper mathematical limits for different hinge arrangements 
in a hinged rigid— beam model are imposed on the elastic— beam equations with 
rigid-body mode shapes, the correct hinged, rigid-beam equations, for each 
particular hinge arrangement, are obtained. 

10. Any three angles can be used to describe the orientation of any 
hinged, rigid body as long as the constraints that relate the rotations about 
the hinges to the variables in the derivation are properly accounted for. Of 
course, it may be more natural to use the physical angles as variables. It is, 
therefore, not the transformation sequence that influences the dynamic 
behavior — it is the physical arrangement of the hinges. 

11. The bending curvatures, torsion, and [T] matrix are developed in the 

appendix for six possible sequences of the angles 6, and 0. Since values 

of the torsion and elements of the [T] matrix are fixed at a given instant in 
time, a change of variable may be found from either quantity to relate the 
geometry from any one sequence to that of any other. 
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APPENDIX A 


KINEMATICS FOR OTHER SEQUENCES 


In the text, the development was restricted to only two of the six pos- 
sible sequences of rotations In the appendix we will present expres- 

sions for the [T] matrix, the bending curvatures, and the torsion for all six 
possible sequences of rotations: ( 1 ) 6 ^, 0 ^; ( 2 ) d^; ( 3 ) Cj, 63 , 

^3’ ^4’ ® 4 » ^ 4 » ® 5 > ? 5 > P 5 > ( 6 ) 9g, B 0 j The necessary changes 

of variable to demonstrate equivalence can be formulated from the torsion and 
then the elements of [T] can be shown to be equivalent for all six sequences. 
Alternatively, the change of variable can be formulated from the T 2 3 element 
of [T] and then the torsion and bending curvatures can be shown to be equiva- 
lent for all six sequences. It is the second approach that is employed here. 

We will make use of equations (5) extensively, which we repeat here for 
convenience 


(x + u) = T 33 ^ 

12 

'13 / 


v"*" = T, 


w+ = T, 


(Al) 


These relationships hold for all possible sequences and types of rotations. 
The [T] matrix, the bending curvatures, and the torsion are now presented for 
each of the six sequences. 

( 1 ) matrix is 


[T] = 


c 

61 


“C„ s 


S„ S„ c 




C_ - s_ s„ s. 


1 ^1 9l 3i 01 Cl Cl Bi 01 Bi'^ei 

■''ei^Bi^Ci ®01^1 ~®0i^Ci ■ ^Ci^ei^ei '^Bi'^ei 


The bending slopes are, from equation (Al) 


(A2) 


v+ = c s ; 
61 Cl 


w"*" = s 


Bl 


v' = Cl + 0 (e 3 ) 


w' = 61 + 0(e3) 


(A3) 


The rotation rate vector is 


< = 


^ ~ ■ ®Ci^) ^ + <jj + <1 


(A4) 
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Thus 


K . = 


9^ + = 0 ! + v'*w’ + O(e^) 

Pi 1 


K . 

J 


-etc„ + ^■J^c s = -w" + e^v" + 0(e3)| 


' 1 ' 


1 


1 


(A5) 


K, = 




(2) ^ 2 , ^ 2 .> ®2‘~ matrix is 


[T] = 


c c 

?2 32 


^32^02 


^.2^3, 




(A6) 


-So + s„ s^ c„ -s„ c_ s„ s„ s_ + c„ c 


32 @2 09 c, 3, 0, ?, 3o 00 ^ 


2 ^2 ^2 


2 ^2 ^2 2 ^2 


^2 ®2 


The bending slopes are, from equations (Al) 


v“^ = s 


W+ = c s ; 

?2 32 


= ?o + O(e^) 


w' = 3^ + 0(e3) 


(A7) 


The rotation rate vector is 


K = -3+J -i ct(cg^K - 


1 + 02^- = Ki + <j+<k 
/ i j k 


(A8) 


Thus 


K = 6, - &'ts = 6' - v’w" + 0(e3) 

1 ^ /LQ Z 


" ■^ 2 ^ 02''?2 ^ 
= V62 ^ ^2®02''?2 ^ 


(A9) 
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(3) Qgj ^ 3 *“ matrix is 


[T] = 


c_ + s„ s„ s 


'3 S 


'3 "3 ‘’3 


-c^ s. 


3 ^3 


®r " ®R ®f) ®i 

^3 ^3 ^3 ®3 ^3 ®3 ' 


c c 
®3^3 


“C s + c s s -c s c - s s c c 

^3 °3 ^3 ^3 ^3 ^3 ®3 ^3 ®3 ^3 ^3 


The bending slopes are, from equation (Al) 


3 S ^3 ’ 


v^"=c s -s„s„c 
^3 ^3 


w"^ = s 


^3 63 


v' = C3 - ^363 + 0 (e3) 


w’ = 63 + 0(e3) 


The rotation rate vector is 


Thus 


« . + 0 + (c^^I + s J) - 6+J - K.l + K.J + 


' 5tc„ s. + etc„ - 6' + v"u' + 0(e3) \ 


■36363 303 -3 


"3 ■ + 4=63 ' -“" »3''" + 0("’) 

“k = «3”" + 


(4) 3,^, 0^, The [T] matrix is 


[T] = 


c c - s s s 
^4 ^4 ^4 ®4 ^4 


^4 ^4 ^4 ^4 ®4 


■"94'^ 


^4 4 


^4 ^4 ^4 ®4 ^ 


^4 ®4 ^4 ^4 ^4 ®4 ^ 


-S, 


c c 

64 04 


The bending slopes are, from (Al) 


v"*” = s_ c. 


w"*" = c, s„ + s_ s„ c. 


■4 -4 ?4 04 ^ ’ 


v' = + 0 (e3) 


w’ = 6^ + c^0^ + 0 (e4 


3 

(AlO) 

3 


(All) 


(A12) 


(A13) 


(A14) 

4 


(A15) 
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The rotation rate vector is 




^1 ■ <% - 


- v’w" + 0(e3) 


K. = -Btc c - 0 ^S = -w" + 64 V" + 0(e3) 

J ^ ?4 °4 ^4 


= ct + S, S = v" + 01^ w" + O(e^) 


'k "4 


(5) 65 , C 5 , B 5 .- The [T] matrix is 


- \'‘h Wo. " 


[T] = -s, 


^5 ®5 


^5 ®5 




The bending slopes are, from equations (Al) 


v+ = c s c - Sg s - C 5 - 3565 + 0(e3) 
^5 ^5 ^5 ^5 5 


j+ = c„ s^ s + s c 

^5 ^5 ®5 ^ ®5 


S 5 + ^565 + o(e^) 


The rotation rate vector is 


* ^^(>= 63 '= - 


,i + K . j + Ki It 


etc, + ct«« = el + + 0^"'^ 

5 $5 C5 5 iig 5 

_«■'■- e‘''s = -w" + 0 ^v" + 0 (e 3 ) 

5 5 5 


.+c - e+c s, = v" + eV' + 0(e3) 

O h 0 i,h 0 
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( 6 ) 6 g, gg, ^ 6 '~ matrix is 


[T] = 


W\W 

-\\ \\^\\\ 


-S , 






The bending slopes are, from equation (Al) 




+ 0(^3) 

^ bo 


= ^6 + ^ 6®6 + 


The rotation rate vector is 


K = 0+1 - B+(Cg J + Sg 

b b 


f«6'' ■" ' "d + + “k*- 


(A22) 


(A23) 


(A24) 


Thus 


"l ■ = ®6 - 

'j ■ <=6 - -“" + D,v" + 0(e3) 

b b 6 

^k = 4 - + 06”" + 


I 


(A25) 


Demonstration of Equivalence 

It can be shown that changes of variables constructed from equating three 
elements of the [T] matrix for any two sequences will lead to identical expres- 
sions for bending curvatures and the torsion. Only one example is presented 
for illustration — the equivalence of (1) and (2). For convenience we choose 
to obtain the change of variable from T 23 . From equations (Al) , (A2) , 
and (A 6 ) , the following relations hold: 




33 = <=6.=^ = - \%\ 


1 


(A26) 
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Differentiation of equations (A26) yields 




■ ‘’I®Bi®0j ■ ® 2 (= 02 ‘=e 2 '^“ 82 '‘^ 2 '*e 


" '82^2=82) ■ 


^1 - c? 



c 
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Differentiation of equations (A26) yields 






1 




3ts„ s, =eVc 


1 3 , 0 ^ 


^2 '“'2 


©2 ^2 P 2 > 


Thus 





s^ s 




Substitution of equations (A28) into (A27) yields 


, ^ (,2 ^2 ^2 '^2 
P-, - 




1 - 

S ^2 


,+ /g^c^ c„ - ?ts. S, 


+ S 2 P 2 '’2 

^ ~ I- 

1 


2 ^2 B 


2 P2 


1- c2 s2 
?2 ^2 


, , ^2 ^2 ‘=2 '-2 
01 - — 


1 - C^ 

^2 ^2 


Substitution of equations (A27) through (A29) into equations (A5) yields for 


K . 

1 


, “6,'C 

Ki - 6, 


2 '’2 ^2 ^2 


1 - c2 s? 


C 


2 P2 


4 ®R ®r 

P2 ?2 


ro+ + 

^2 '^2 ^2 ^2 
2 o2 


1 - s 


>2 ^2 


= 0^ - s 3^ 
2 ?2 ^ 


(A29) 


which is identical to the expression for Kj^ in equations (A9) . The bending 
curvatures also transform identically. Similarly, the other four 
sequences (3)-(6) can be shown to be equivalent to (1). The changes of 
variable to second order may be determined by inspection, which proves that 
all six sequences are also equivalent to second order. 
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